Current-carrying cosmic string loops 3D simulation: 
towards a reduction of the vorton excess problem. 
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The dynamical evolution of superconducting cosmic string loops with specific equations of state 
describing timelike and spacelike currents is studied numerically. This analysis extends previous 
work in two directions: first it shows results coming from a fully three dimensional simulation 
(as opposed to the two dimensional case already studied), and it now includes fermionic as well as 
bosonic currents. We confirm that in the case of bosonic currents, shocks are formed in the magnetic 
regime and kinks in the electric regime. For a loop endowed with a fermionic current with zero-mode 
carriers, we show that only kinks form along the string worldsheet, therefore making these loops 
slightly more stable against charge carrier radiation, the likely outcome of either shocks or kinks. 
All these combined effects tend to reduce the number density of stable loops and contribute to ease 
the vorton excess problem. As a bonus, these effects also may provide new ways of producing high 
energy cosmic rays. 

PACS numbers: 98.80.Cq, 11.27.+d 



I. INTRODUCTION 

Superconducting cosmic strings are topological de- 
fects U that may be formed at a phase transition in 
the early universe Contrary to their usual coun- 

terpart known as Goto-Nambu strings Q, these possess 
a nontrivial internal structure because of conserved cur- 
rents that can flow along them S 0, M. The dynamics 
of superconducting cosmic strings is, because of the cur- 
rent structure that has to be taken into account, much 
more complicated. In order to describe their evolution, 
the so-called elastic string formalism, thanks to which 
the micro-structure can be integrated over to yield the 
macroscopic behavior, was set up by Carter j8|. The only 
requirement is the knowledge of the equation of state re- 
lating the energy per unit length to the string tension. 
This is equivalent to giving a single Lagrangian function 
C(w), depending on a so-called state parameter w, out 
of which the dynamical equations can be derived. 

Various macroscopic Lagrangians for describing super- 
conducting cosmic strings have been proposed ||, that 
correspond to different assumptions about the internal 
structure. In principle, the current along a cosmic string 
exists because some extra degrees of freedoms (e.g. parti- 
cles are trapped in the defect core. Up to now, there 
are essentially three models that have been developed to 
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yield analytic equations of state. Those can be called 
respectively transonic, bosonic and zero mode fermionic 
carrier models. 

The transonic model describes for instance wiggly 
Goto-Nambu strings Q. The trapped degrees of free- 
dom in this case are not particles, but rather transverse 
string excitations that play the same role. It is self-dual, 
meaning that the equation of state algebraic form does 
not depend on whether the current is timelike or space- 
like, and it is transonic because the velocities of trans- 
verse and longitudinal perturbations happen to be equal. 
Among its numerous advantages, one finds that the tran- 
sonic model can be exactly solved in the case of flat space, 
and the string loop motion in this model is always sta- 
ble. It also describes the motion of an ordinary string 
in a five-dimensional space-time, a la Kaluza-Klein [ pd| . 
In this last case, the state parameter represents the fifth 
coordinate of the string location, seen by projection in 
the four dimensional space-time. 

Superconducting cosmic string may also involve ac- 
tual particles, the simplest case being that of fermions 
trapped in the defect core in the form of zero modes, as 
originally proposed in 1985 by Witten ^ to introduce su- 
perconductivity. Under very general circumstances, the 
fermion condensate can be shown Jl2| to rapidly reach 
a zero temperature distribution, so that the integrated 
model does not depend on anything else but the string 
characteristic mass scale, m say. The relevant equation 
of state in this case is then again self-dual, with the en- 
ergy per unit length and tension adding up to a constant 
at all times. 

Bosonic current-carrying cosmic strings are described 
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by an equation of state which is not self-dual. Such 
models require two different masses the string scale 
m and a mass characterizing the current, say, usu- 
ally taken to be that of the carrier itself. Such mod- 
els have been examined with full details || fjj. Even 
though a bosonic condensate can be treated as a sin- 
gle classical field, contrary to it fermionic counterpart, 
the corresponding strings look rather more involved than 
the fermionic current-carrying ones. Indeed, it turned 
out that two equations of state are required to describe 
bosonic current-carrying cosmic strings, one for each of 
the available regime, namely spacelike or timelike. What 
previous analysis revealed is that even in this case, quan- 
tum effects should be taken into account, due to various 
possible instabilities that could be triggered by the loops 
evolution fL3| , 

There are results on the evolution of superconducting 
string loops in one and two dimensions |D| |l^| for those 
equations of state corresponding to strings with bosonic 
currents. These works have shown the appearance of var- 
ious singular behaviors as the dynamical evolution can 
drive parts of a given loop outside of the domain of va- 
lidity of the elastic string description. For instance, the 
loop can develop shocks |l2 , 15 or fold on itself in compli- 
cated shapes. The first question that could then be asked 
regarding these effects is how much of these is due to the 
restriction to one or two dimensions ? In other words, 
are these results real or do they merely arise because of 
some projection effect. In the present paper, we have ac- 
cordingly extended previous studies to the simulation of 
a superconducting cosmic string loop with bosonic cur- 
rents in three dimensions. We find that all the singular 
behaviors observed in two dimensions generalize to three 
and cannot therefore be assigned to some projection ar- 
tifact. 

Furthermore, we also investigated the behavior of a 
loop with fermionic currents, using what was recently 
shown to be the relevant equation of state for such 
strings (l^, i.e. the fixed trace model. In this case, 
although dynamical evolution can also drive the loop out- 
side of the domain of validity of the elastic string de- 
scription and fold it in complicated shapes, shocks never 
develop. All these singular behaviors, both for fermionic 
and bosonic carriers, are interesting because, they are ex- 
pected to lead to charge carrier emission, a mechanism 
that has many cosmological consequences. 

The dynamics of the string loops under consideration 
in this paper demands to be applied to a network |l7| of 
such strings in order to decide on important issues such as 
vorton formation |H| and evolution jf{| . Another related 
possibility concerns the up-to-now mysterious highest en- 
ergy cosmic ray events |20 . These could be explained 



in terms of topological defects such as superconducting 
cosmic strings, and indeed there have been many such 
proposals |2l| . In the case at hand, our conclusions seem 
to imply that much more particles are released during 
the normal evolution of a current-carrying cosmic string 
loop than during that of an ordinary loop. Assuming 



these particles to have enough energy to contribute to 
the highest energy cosmic ray, this paves the way to a 
new class of models that may not suffer from the nor- 
malization constraint B^j . 

The article is articulated as follows. In the next section 
we derive the equations of motion in three dimensions, 
and express them in an appropriate way for efficient nu- 
merical resolution. In section III, we present the equa- 
tions of state for superconducting cosmic strings with 
bosonic and fermionic currents, while the results of the 
simulation are exhibited in section IV. For a loop with a 
bosonic current, we exemplify cases where the loop de- 
velops shocks, kinks or folds on itself in much the same 
way as in the two dimensional case. The analysis is then 
extended to superconducting loops with fermionic zero- 
mode current-carriers. In this case, we find that in some 
instances, the loop shows regions of discontinuous cur- 
vature (kinks) but we do not find shock waves as in the 
bosonic situation. We also observe that in most of the 
cases we have investigated, the dynamical evolution of 
a loop, whether with bosonic or fermionic current, will 
drive it out of the elastic regime. This is yet another 
confirmation that the quantum (microscopic) effects that 
are taking place in a cosmic string loop are almost always 
not negligible, and even more so in the cosmological set- 
ting for which those effects might modify drastically the 
model predictions |^|. We discuss this particular point 
in the concluding section. 



II. EQUATIONS OF MOTION 

For a current-carrying cosmic string, the equations of 
motion can be expressed as the conservation of the stress 
energy tensor and the equation of state. The stress en- 
ergy tensor of a cosmic string can be expressed in diago- 
nal form as H : 



(1) 



where and v M are the two orthogonal, respectively 
timelike and spacelike unit eigenvectors 1 , which describe 
the string worldsheet, and U and T are the two cor- 
responding eigenvalues identified respectively with the 
energy per unit length and tension of the string. The 
equations of motion stem from the conservation of the 
stress energy tensor, and they can be split into two pairs 
of intrinsic and extrinsic dynamical equations. The in- 
trinsic dynamical equations are obtained by projection 
along the string worldsheet while the extrinsic dynam- 
ical equations are obtained by projection perpendicular 
to the string worldsheet. The intrinsic equations reduce 
to two current conservation laws, one timelikc 



Vp(i/u p ) - 0, 



(2) 



-1, u u vV 
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and the other spacelike 



V p (Ai«") - 0, 



(3) 



where V p = 77p <T V CT is the covariant derivative on the 
string worldsheet, obtained by projecting the usual co- 
variant derivative on the worldsheet with the projection 
operator r\ p a = —u p u a + v p v a , and v and [i are defined 
through the relations 



[IV 



U-T, 
dU 



(4) 
(5) 



and can be understood, respectively, as a number den- 
sity variable and its associated effective mass variable or 
chemical potential. 

The intrinsic equations can be equivalently re- 
expressed as irrotationality equations, respectively as 



and 



e ap V a (fj,u p ) = 0, 



(6) 



(7) 



where e ap is the antisymmetric tangent element tensor of 
the string worldsheet, given by 



(8) 



in terms of the tangent vectors u and v. 

For a superconducting cosmic string, the conserved 
numbers associated with the two conserved currents cor- 
respond to the charged current trapped in the string and 
the winding number of the string. When a loop is con- 
sidered, these translate into two integer-valued charges, 
N and Z, say, once integration around the loop is per- 
formed. 

The extrinsic equations of motion can be written as 



L p p {Uu u V u u p - Tv v V v v p ) = 0, 



where _L M „, defined by 



- 1 - v — y v '/ v — y u 



(9) 



(10) 



is the orthogonal projection operator. 

The necessary quantities for describing a cosmic string 
loop are completed by the knowledge of two other im- 
portant parameters, namely the speeds of transverse and 
longitudinal perturbations along the string, written re- 
spectively as c T and c L , and given by 



T 



dT 
dU' 



(11) 



These velocities are useful to explore the stability of a 
loop at equilibrium |p6[ . Furthermore, they are helpful 
since they allow us to define what we call elastic regime: 
the string dynamics is in the elastic regime if transverse 
and longitudinal perturbations are stable. This indeed 



ensures that the dynamical model does not break down. 
The velocities (|ll| ) should also be asked to both be less 
than unity in order to avoid superluminal propagation, 
but this, in practice, never actually happens so that it is 
not really a constraint. To summarize, we require that, 
along the loop trajectories, we have 



0<< <1, 



(12) 



constraints that are not always trivially satisfied. 

We now have to fix the gauge to determine unambigu- 
ously the relevant degrees of freedom. The dynamical 
equations are solved numerically choosing as unknowns 
the string worldsheet coordinates a; M (i,-0) as in and 
as parameters t, the time coordinate, and ip, a spacelike 
curvilinear coordinate chosen as the potential associated 
with the irrotationality equation (||). With this choice, 
u p , v p and v are expressed as: 



,p - 



■ v i> 



c 

n 



,((3x p + CV"), 



(13) 
(14) 

(15) 



with dots and primes respectively denoting derivatives 
with respect to t and ip and we have adopted the following 
definitions 



C = y/-±Px p , 

[3 — x x p - 



(16) 



(3 2 + C' 2 C 2 



The equation of motion (^|) is automatically satisfied by 
the choice of its associated potential ip as variable, and 
the remaining equations of motion (Q) and ^ become 
respectively 

v • [(n 2 - /3 2 c 2 )x - C^x" - 2/?C 2 c 2 x'] = 0, (17) 

where v is the spatial part of i> p , and 

wxi,2 • l(n 2 - P 2 c 2 T )x - C 4 c'x" - 2/3C 2 c 2 x'] - 0, (18) 

with wn^ the spatial parts of two independent 
quadrivectors orthogonal to the worldsheet, which we are 
choosing as 



1 



W±l 



V z - yz 



— (y'z-yz',0,-z',y'), 



(19) 



and 



W± 2 



y z-yz 



7(0, y'z - yz', xz' - zx', yx' - xy'). (20) 
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Then, the equations of motion can be solved to find that 

x(t,%l>),y(t,ip) and z(t,ip) satisfy 



s{z'F + K), 



(23) 



x = \{x'F + G), 
v = \{y'F + l), 



(21) 

(22) where the functions F, G, I and K are given by: 
I 



2^,2 



F 



G = 



I = 



K 



n 2 — (3 2 c 2 

n 2 — P 2 ^ 

C 2 c 2 
n 2 — /? 2 c 2 

n 2 — f3 2 c 2 



C 2 (cix" + c 2 y" + c 3 z") + 20(c L x' + c 2 y' + c 3 z') 



(x'y" - x"y')c 2 + (x'z" - x" ' z')c 3 

{y"z' - z"y')c 3 + (y"x' - x"y')ci 

(y"z' + z"y')c 2 + (z"x' - z'x")ci 

I 



+ 2/3 
+ 2/3 
+ 2/3 



{x'y' ~ x'y')c 2 + {x'z' - x'z')c 3 
{y'z' - y'z')c 3 + {y'x' - x'y')c x 
{y'z' + z'y')c 2 + (iV - z'x')c\ 



(24) 
(25) 
(26) 
(27) 



with 

ci = x' + yA 3 + zA 2 , 

c 2 = y'-xAx+zAa, (28) 

c 3 = z' - xA 2 - yA 3 , 

and 

Ai = y'x-yx', 

A 2 = z x — zx , (29) 
A 3 = z'y - zy'. 

As a test of the deviations of the numerical approxima- 
tion from the exact solution it will be convenient to calcu- 
late the (conserved) total energy and angular momentum 
of the system given, with our choice of unknowns, by {I 
being the string arc length) Jl3| 

E T = jdlT™, (30) 

which reads, explicitly 

E T = I d^fll(n 2 - (31) 

and the angular momentum Q 

J T =ZN = * fd^^{l-c 2 T ), (32) 

where Z is the particle quantum number, N is the wind- 
ing number and ^ is the parametric length of the loop in 
the ip coordinate (which is conserved, actually being one 
oiNoiZ). 



We shall assume in what follows that the dynamics 
obtained by means of the numerical solution is accurate 
as long as both these quantities are conserved at least at 
the few % level (worst case during the development of a 
shock) all along the trajectory. 

III. EQUATION OF STATE 

The equation of state describes the substructure of the 
current-carrying cosmic string considered by relating its 
internal parameters, for instance its energy density and 
tension. There are three interesting current-carrying cos- 
mic string models: the transonic model and the bosonic 
and fermionic superconducting models. The transonic 
model is used to describe the macroscopic dynamical be- 
havior of a wiggly "Goto-Nambu" cosmic string and car- 
ries a current associated with the wiggles jl0| : it permits 
to neglect the dynamics of the strings on small scales by 
integrating over these scales, effectively raising the stress- 
energy tensor degeneracy, and thereby to introduce an 
effective current. This model also possesses the nice fea- 
ture of being algebraically solvable in flat space , and 
thus can provide useful analytical examples with which 
one can compare the accuracy of a simulation code. The 
other models we shall deal with are used to describe cos- 
mic strings in which a current of particles (either bosonic 
or fermionic) has condensed. 

For superconducting strings with a bosonic current, 
the model requires two mass scales, namely the cosmic 
string mass scale m, which is the string forming symme- 
try breaking energy scale, and to* which is the mass scale 
of the current carrier. In this case, one can write down 
two different equations of state depending on whether the 
conserved current along the string is timelike or space- 
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like |(| . These two regimes of the model are called respec- 
tively electric and magnetic in accordance with the situ- 
ation where the current is coupled to an electromagnetic 
field 0. Indeed, in this CM SC. c.l timclike current would 
induce an essentially electric field (in the sense that, as 
E 2 > B 2 , the magnetic component can be framed away), 
while a spacelike current leads to a magnetic-like field. 
Using a toy model which produces such strings ^| , the 
equation of state for U as a function of T in the magnetic 
regime is M 



2 m Z m * 
U = m 2 + — - - — 



m% 
Iff 



T, 



(33) 



and the transversal and longitudinal speeds are given by 
1 2 m 2 (\ + v 2 ) 2 -v 2 {\-v 2 ) 



c 2 = 



1 + ^ 2 
1 - 3^ 2 

l + vi : 



2m 2 (1 



(34) 
(35) 



where i/* = v/m*. In this case, c 2 can become negative 
so that the string may become unstable with respect to 
longitudinal perturbations and thus leave the domain of 
elasticity. In the electric regime, the toy model yields an 
equation of state of the form 



U = T 



-T)/r, 



1]. 



and the two perturbations velocities are 

2 _ [2(m 2 /m 2 ) + ln(l - X*)] (1 - X*) 



[2(to 2 /to 2 ) + ln(l - X,)] (1 - X,) + 2X* ' 



where 



and 



X, 



2v 2 



\-y/E 



21/2 



c 2 = 



v/1 + 4^ 2 



(36) 



(37) 



(38) 



(39) 



In this regime, it is clear that only c 2 may become nega- 
tive, at which point the string may become unstable with 
respect to transverse perturbations. 

The final model of current-carrying cosmic strings 
which we shall consider here concerns the case in which 
fcrmions may be bound to cosmic strings in the form of 
two lightlike currents of zero modes propagating in op- 
posite directions. The resulting total fermionic current 
can then be timelike, spacelike or lightlike. Because the 
fcrmion condensate tends to a zero temperature distri- 
bution, the equation of state is self-dual, that is of the 
same form whether the fermionic current is timelike or 
spacelike, given by |l2|, |l6f 



U + T = 2m 



(40) 



which gives 
U 



2 ' 



T 



v 

y 



(41) 



where m is, as before, the string scale. Here, the pertur- 
bation velocities are simply 



2m 2 2m 2 — v 2 

- 1 = 



U 



2m 2 



1. 



(42) 
(43) 



In this case, the string can leave the domain of elasticity 
only if the velocity of transverse perturbations c 2 can 
become negative, i.e. if \u\ > \[2m. We shall see that 
this indeed happens. 



IV. EVOLUTION 

We studied the evolution of loops in three dimensions 
after having reproduced, as a check, exactly the two- 
dimensional results of Ref. |u|. Then, the evolution of 
loops in three dimensions was investigated with the var- 
ious superconducting equations of state discussed above. 
The class of initial configurations with which we have 
been concerned was obtained by perturbing the circular 
equilibrium solution in the following way: 



x(t = 0,-0) = i?cos 
±(t = 0,ip) = 



/ 4, 



\Rr/vo 



c T . ( ip 
sm f 



y(t = o,v) 

z(t = 0,ip) 
z(t = 0,0) 



1 + e \ Rjv 



Re sin 



R'yvo 



(■-,. I Ip 
■ cos 



1 + e 



Riv 



= At 



R~/v 



—B sm 

V R-yvo J 



(44) 
(45) 
(46) 
(47) 
(48) 
(49) 



where vq and 7 = (1 — x 2 — y 2 — z 2 ) 1 / 2 are respectively 
the current and the Lorentz factor of the circular equi- 
librium solution, e is the ellipticity of the configuration, 
ijj is the spacelike coordinate, chosen to vary from zero 
to 2irR r yvi) using the gauge freedom in the worldsheet, 
e is a dimensionless parameter which measures the de- 
viation in velocities from the circular equilibrium state 
obtained when e = 1 and e = 0, and i is an integer phase 
parameter for the transverse z motion. Note that the 
amplitude A of the perturbation in the transverse (with 
respect to the string loop plane) direction, and that of the 
associated velocity B are independent since those pertur- 
bations have been shown pq] to be absolutely stable. 

For a given equation of state, the stability of the per- 
turbed circular equilibrium solution can easily be de- 
rived |2q]. The primary purpose of this work is to find 
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out the fate of the unstable configurations. Before study- 
ing such configurations within the framework of bosonic 
and fermionic superconducting equations of state, we first 
checked the program against stable ones with the tran- 
sonic equation of state. In this case, everything runs 
smoothly and there is nothing particular to note about 
the loops dynamics. Thus, no evolution of transonic loops 
is shown (see however Ref. j[3| for a discussion of these 
configurations; the three dimensional simulation now pre- 
sented merely confirms this reference results as we were 
able to check that they were not artifacts of the two di- 
mensional projection). 

More interesting is the case of the bosonic current 
where singular dynamics where observed in two dimen- 
sions The analysis in this case is presented next. 
It is then followed by an investigation of the fermionic 
current case which is completely new. 



A. Bosonic Currents 

For superconducting strings with bosonic currents, we 
confirm that all the singular behaviors observed in two 
dimensions, i.e., shocks in the magnetic regime, kinks in 
the electric regime, folding loops and exit from the elastic 
regime in general, persist in three dimensions ; they are 
not projection artifacts and represent the real physical 
time evolution of superconducting cosmic string loops. 
In the following, we show examples that illustrate each 
case in turn. 

A shock wave may be formed at points along the string 
worldsheet, and according to previous results H H we 
further confirm that a shock can only take place in the 
magnetic regime. An example of such a shock in three di- 
mensions, using equation of state (|3^) is shown in Fig. |l| 
on which is also represented the space variations of the 
state parameter = v/m for different instants in time 
for a transversely perturbed loop with I — 3 ; note that 
such a loop, from the transverse perturbation point of 
view, should be strictly stable. In this case, the geomet- 
ric evolution does not reveal any visible change in the 
loop, but the time evolution of the state parameter vl as 
a function of the tp parameter shows the development of 
shock fronts along the string worldsheet. At these points 
the current in the string becomes discontinuous and the 
code cannot handle the evolution anymore, as shown for 
instance by the fact that the total energy is no longer 
conserved. Although such transverse modes in principle 
would have been expected to be stable, it turns out, with 
the values of the parameters chosen, that nonlinear ef- 
fects rapidly become dominant, effectively coupling the 
transverse mode to the longitudinal one. That this really 
occurs is made clear by the fact that the number of shock 
fronts generated during this evolution is 2£, thus directly 
related to the phase parameter £. 



Another important effect that may take place along a 
superconducting loop is the appearance of kinks. These 
are regions where the curvature of the string becomes 
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FIG. 1: Time evolution of the state parameter v, along the 
string worldsheet for an elliptic configuration of a loop with 
a bosonic current in the magnetic regime, having large ellip- 
ticity e = 0.6 and velocity parameter e = 1, and with a phase 
parameter £ — 3 (and A — 3, B — —0.2). In the last panel 
the variation of the energy shows two important things. On 
the one hand the small peaks in the graph are an evidence 
that the loop size decreases because the initial velocity has 
been reduced. On the other hand, nearly at the end of the 
simulation, the rare variation of the energy gives a signal of 
the shocks produced along the loop string. The relative varia- 
tion of the total energy, shown on the last panel, is still rather 
small even as the simulation reaches its very end. The situa- 
tion is similar for the angular momentum, and the simulation 
eventually breaks down. 



discontinuous. We found that giving an extra dimension 
for the loop to evolve into did not suppress the appear- 
ance of kinks. Fig. || shows the evolution of a loop in 
the electric regime-and therefore using the equation of 
state given by Eq. (gfj). In this case again, the simulation 
shows that 2£ kinks are being formed at points along the 
loop. To illustrate this point further, we plotted in Fig. || 
the curvature of the loop as a function of the spacelike 
string internal coordinate tp. This is, in general, given 

by II 



i.e., in the case at hand, 

K 2 = (w p xl K p ) 2 + (w^ 2 K p f, 

where 



(50) 



(51) 
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(n 2 - /3 2 ) (yz' - zy') + 2/3C 2 {y'z> - z'y>) + C 4 (yV' - z'y") 



(52) 



(n 2 - /3 2 )(A 3 x - A 2 y + A x z) - 2/?C 2 (A 3 i' - A 2 y> + A 1 z') - C 4 (A 3 a;" - A 2 y" + A lZ ") .(53) 




FIG. 2: Configuration evolution in the electric case with 
bosonic currents with the same ellipticity as in Fig. [j], e = 0.1 
and A = 3 (B = 0). One clearly sees the development of six 
kinks corresponding to twice the phase parameter £ — 3. In 
this case, the peaks seen in the energy are formed at the same 
time as the kinks are present. 



It is clear in the figure that where kinks appear, large 
variations of the curvature are developed, showing that 
the kinks are not gauge artifacts but really occur. As 
expected, the location of the kinks corresponds to the 
points with large curvature. In general the effects found 
in the two dimensional planar case are still present in 
three dimensions. They are however usually affected (or 
even initiated) by the transverse perturbation in z. 

In some cases, the string loop shows points where in- 
tercommutation, which is not taken into account in our 
code, should be occurring. A systematic analysis shows 
that a loop in three dimensions tends to fold on itself in 
any regime, magnetic or electric. This is similar to what 
happens to Goto-Nambu loops for which it was even ar- 
gued that this intercommuting process could be the lead- 
ing process for a string network to loose energy |2^, p5[ . 
The conclusions that can be drawn from this fact for a 
superconducting string network are therefore the same, 
and will be discussed later. In Fig. |] is shown an exam- 
ple of such a missed intercommutation for a loop in the 
electric case. 

Another possibility that is open for a loop to leave the 
region where the elastic description is valid is to locally 
have one of its squared perturbation velocity, c 2 in the 



r 




2 4 6 8 



FIG. 3: Invariant magnitude K = V K 2 of the curvature, from 
Eq. (|5l]), for the configuration shown in Fig. |ij 

magnetic regime, or c 2 in the electric regime, becoming 
negative. In the latter case, the tension reaches nega- 
tive values and the string actually behaves like a spring. 
In Fig. H is shown the evolution of the longitudinal per- 
turbation velocity as a function of the spacclikc parame- 
ter along the string loop at different instants of time in 
a typical example for which the dynamical evolution of 
the loop drives it out of the elastic regime when c 2 be- 
comes less than zero. Whatever happens past this time is 
presumably not accounted for correctly by the program 
which, being an effective macroscopic description, does 
not recognize the instability at the microscopic level. 

All these results are reminiscent of what had previ- 
ously been obtained in the two dimensional case |l3[ . 
There now is no more ambiguity about them being pro- 
jection artifacts. The conclusions derived in the 2D case 
therefore still apply. 

B. Fermionic Current 

In this section, we consider current -carrying cosmic 
strings with fermionic currents, obeying the equation of 
state (pp|). The behavior in this case follows somewhat 
that observed in the bosonic case: the loop can fold on 
itself and leave the elastic domain. We shall exemplify 
both cases in turn, but first let us remark the important 
following point. Since in the fermionic case the equation 
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FIG. 4: A typical string loop configuration showing points 
where intercommutation should take place, with 1 = 2, el- 
lipticity e = 0.3, A = 1, B = 0.1, and vt = 0.1. At these 
points, two loops should be formed and the subsequent evo- 
lution given by the code is no longer physically relevant. 



of state is self-dual, there cannot be any qualitative dif- 
ference between the magnetic and the electric regimes at 
the level of the equations of motion. However, the two 
perturbation velocities are not equal and the string may 
become unstable only with respect to transversal per- 
turbations. As a result, the dynamical evolution of the 
string can lead it to develop only kinks 2 , although in this 
case we cannot say that the effect is exclusive of the elec- 
tric or magnetic regime. In Figs. |^ and |t], we show the 
time evolution and the curvature (snapshots at different 
times) for a loop with fcrmionic current, again with phase 
parameter t = 3 and initial value of the state parameter 
set as v^q =0.1. The appearance of kinks is seen as the 
curvature of the string becomes discontinuous. When the 
kink is formed the total energy is no longer conserved and 
the simulation ends. 

Similarly to the case of cosmic string loops endowed 
with bosonic currents, loops with fermionic currents 
can fold on themselves, exhibiting points where inter- 
commutation should take place, and with very much 
the same rate of occurrence in the evolution. Self- 
intercommutation for loops loosing their energy again 
seems to be in this case a favored mode. 

The string loop can also leave the domain of validity of 
the elastic formalism when c 2 becomes negative. Fig. ^ 
depicts the time evolution of the transverse velocity c 2 
up to the point where it vanishes. The velocity here is 



2 Note that since shocks can only form in the magnetic regime and 
as in this self-dual situation, the electric and magnetic regime are 
undistinguishable, shocks cannot form in the fermionic current 
case. 
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FIG. 5: Time evolution of the squared perturbation veloc- 
ity, c 2 , in the magnetic regime for a bosonic superconducting 
loop. The dynamics drives the loop out of the elastic regime, 
as seen on the last evolution panel, in which one finds a re- 
gion having c 2 < 0. The last panel shows the total energy 
integrated along the loop as a function of time: its constancy 
up to the point where becomes negative reveals that noth- 
ing particular took place and the physical description can be 
trusted. Past this point, the energy is very badly conserved 
and the program, as expected, has difficulties handling the 
problem. This figure is for A = 2 and B — 0. 



shown as a function of the spacelike string parameter ip. 
The same remarks as for bosonic strings apply, in partic- 
ular concerning the non conservation of the total energy 
whenever the elastic regime is left. Note here the appear- 
ance of a numerical instability that propagates along the 
string and that is seen in the form of (unphysical) oscil- 
lations. 



V. DISCUSSIONS AND CONCLUSIONS 

This work is the direct follow-up of previously dis- 
cussed two dimensional string loop simulations |i~3] ]. It 
extends this old work in two directions. First it now 
provides results about the dynamics of such string loops 
in three dimensions, thereby answering the question as 
to whether the observed effects were actually physical or 
mere projection artifacts. This new simulation clearly 
shows that indeed the effects presented in Ref. |l3j are 
physical. Another extension of Ref. jl3| that is done here 
concerns the equation of state, which can now describe 
fermionic |T^ ] currents as well. 

In the case of strings endowed with spacelike bosonic 
currents, it was shown that most of the loops develop 
shock. These shocks are seen as discontinuities in the 
current and occur because the squared longitudinal per- 
turbation velocity c 2 can become negative. In the electric 
regime of timelike currents, it is the squared transverse 
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FIG. 6: Configuration evolution for fermionic equation of 
state. In this case, the shape of the loop changes until kinks 
are formed along the string, at which point the elastic string 
description ceases to be sufficient. Quasi kinks are formed 
first, that are subsequently rounded off, until the actual kinks 
really form. Here, and in the following figure, the parameters 
A — 3 and B = 0.2 have been used. 
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FIG. 7: Time evolution of curvature against the point labeling 
function ip for the configuration of the Fig. [| The curvature 
in this example varies almost by two orders of magnitude. As 
in the bosonic case, the variations of the energy give a signal 
of the kinks formation. 



perturbation velocity c^, that can vanish and reverse its 
sign. As it turns out, this implies that kinks, namely 
regions of discontinuous curvature, develop. 

Fermionic strings, being described by an equation of 
state for which only the squared transverse perturbation 
velocity can become negative, are subject only to the 
possibility of building kinks, and this in a way that is 
independent of the spacelike or timelike nature of the 
underlying current since the relevant equation of state is 




FIG. 8: Time evolution of the speed of transversal perturba- 
tions for a cosmic string loop with fermionic currents (with 
A — 2 and B = —0.2). The dynamical evolution of the loop 
drives it out or the elastic regime, where < 0. 



self-dual. 

In both cases, fermionic and bosonic currents, and for 
whatever kind of current, the loops tend to fold on them- 
selves, generating contact points where intercommuta- 
tion should occur. The almost systematic occurrence of 
such configuration leads us to conclude that one may 
divide superconducting string loops into essentially two 
species: those loops which tend to divide into smaller 
loops, and the other ones that are identified with the 
proto-vortons of other authors [|l9| . Our results tend to 
indicate that the latter category is much less populated 
than the former, as was already postulated in earlier in- 
vestigations Jigf ] . It also implies that these intercommu- 
tations will be, in the case of superconducting cosmic 
strings, much more frequent than in the case of ordinary 
strings. 

Now when a string loop self-intercommutes, the topo- 
logical stability that confines the Higgs and other parti- 
cles inside the defect is raised and some of these particles 
are expelled away. Although stable while contained in the 
defects, these particles are often unstable. In particular, 
strings are expected to form at the Grand Unified (GUT) 
scale, and the relevant Higgs and gauge fields to decay 
therefore almost instantaneously by cosmological stan- 
dards. These particles will then initiate showers of lower 
mass ones dividing the original energy. This is the mech- 
anism that is at work in most Top-Down scenarios |27j 
of Ultra High Energy Cosmic Rays (UHECR) whose ori- 
gin is still mysterious p0[ . Ordinary cosmic strings do 
not produce much of these particles, because intercom- 
mutations are not that frequent in Goto-Nambu strings 
since most of the energy of the network in that case is ex- 
pected to be lost through gravitational radiation 0]. In 
fact p2| , the resulting flux is found a disappointing ten 
orders of magnitude below the observed one. In the case 
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of current-carrying strings, loops have a higher probabil- 
ity of self-intercommuting, so the mechanism is enhanced 
accordingly. 

The proto-vorton case now takes us to the rare config- 
urations that do not intercommute frequently. Note at 
this point that the scarcity of these states is not enough 
to ensure the vorton excess problem to be altogether alle- 
viated. Indeed, for GUT scale vortons, the ratio between 
proto-vortons and doomed loops should already be much 
less than 10~ 20 in order for the present-day vortons not 
to dominate the Universe |l9j if the resulting vortons are 
stable over cosmologically relevant timescales. 

Gathering all the results obtained over the recent years, 
we now find the following: for a realistic underlying mi- 
croscopic model, the equation of state is such that at 
least one of the perturbation propagation speed can be- 
come imaginary, leading to possible instabilities. The 
main result of the present paper is to show that these 
instabilities will systematically develop in every loop as 
it evolves. As a result, the elastic treatment turns out to 
be valid for limited periods of time and one must resort 
to the full microscopic description in between. 

The elastic description breaks down in two possible 
ways, namely through the formation of kinks or shocks. 
In both cases, the meaning is that the microstructure 
should be taken into account. In particular, close to these 
points, it is no longer appropriate to consider the string 
as effectively two dimensional and we must consider the 
effects that are due to the finite thickness. 

First of all, kinks cannot exist if the string is a tube 
instead of a Dirac line as it would be impossible for in- 
stance to match the Higgs field and its derivatives, as is 
however required by the underlying field theory, at the 
points where the curvature is discontinuous. Therefore, 
the string must somehow find a way of smoothing, which 
means in that context loosing energy. The only non kinc- 
matical (in the sense of necessarily breaking the elas- 
tic description) way of realizing that is through particle 
emission, and we are led to conclude that, as in the case of 
intercommutation, but for completely different reasons, 
some particles (current carriers, but also Higgs and gauge 
fields) will be radiated away as the kinks develop. Note 



also that the current-carrier is indeed likely to get out of 
the string since, as the string curvature increases, the car- 
rier momentum acquires an orthogonal-to-the-worldsheet 
component thanks to which the wavefunction is spread 
over larger distances, thereby effectively increasing the 
probability for the particle to tunnel out. This mech- 
anism can ease the vorton excess problem while in the 
same time providing more efficiency for UHECR produc- 
tion through Top-Down scenarios. 

Shocks develop only for bosonic spacelike currents. In 
this case, the state parameter which would end up 
being discontinuous were the shock actually occurring, 
is well represented by the collective momentum of the 
trapped particles || |?J , this momentum being always less 
than the rest mass of the carrier particles themselves in 
order for the confinement to take place. As we argued 
before Jig ], it is rather dubious that actual shock could 
form in cosmic strings, as the shocks are always accom- 
panied, either followed or preceded, by a transient region 
in which c 2 < 0. However, in this region, the momentum 
of the carrier is higher or at least of the same order as its 
rest mass. Therefore, it is energetically favored to some 
of the trapped particles to move out of the string. 

In all the situations which we have studied, we found 
that the expected dominant effect was particle radia- 
tion. We are therefore led to the conclusion that an 
elastic treatment of cosmic strings is not completely ap- 
propriate. It should be supplemented with some semi- 
classical, non conservative, forces. The dissipative ef- 
fect then introduced can serve to ease the vorton excess 
problem, while providing new interesting ways to produce 
UHECR. 
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